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Quantum  transport  theory  for  semiconductor  superlattices 
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t Paul-Drude-Institut  fur  Festkorperelektronik,  Hausvogteiplatz  5-7, 

10117  Berlin,  Germany 

$ Ioffe  Physico-Technical  Institute,  St  Petersburg,  Russia 

Abstract.  The  Kadanoff-Baym-Keldysh  non-equilibrium  Green’s  function  technique  is  used  to 
study  quantum  transport  within  the  lowest  miniband  of  semiconductor  superlattices  (SLs)  under 
high-electric  fields.  Both  intra-collisional  field  effects  (ICFEs)  and  lifetime  broadening  are  taken 
into  account.  An  inevitable  extension  of  the  generalized  Kadanoff-Baym  ansatz  is  discussed.  The 
collisional  broadening  of  electro-phonon  resonances  due  to  scattering  on  impurities  is  calculated. 
Even  if  the  mean  impurity  scattering  strength  is  considerably  smaller  than  the  miniband  width,  the 
oscillatory  current  anomalies,  which  result  from  ICFEs,  are  completely  smeared  out. 


1 Introduction 


If  the  Bloch  frequency  £2  = eEd/h(E  is  the  electric  field  strength  and  d the  SL  period)  is 
larger  than  some  effective  scattering  rate,  1 /r,  carriers  confined  to  the  lowest  miniband  are 
expected  to  be  Bragg  reflected  before  being  scattered  by  phonons  or  imperfections  in  the 
crystal.  This  gives  rise  to  Bloch  oscillations  and  to  the  formation  of  a Wannier-Stark  (WS) 
ladder.  The  electric  field  induced  WS  localization  results  in  non-analytic  resonant-type 
anomalies  in  the  current-voltage  characteristic  (I— V),  known  as  electro-phonon  resonances 
| ].  Whereas  these  current  anomalies  have  been  observed  in  narrow  band  semiconductors, 
there  are  no  experiments  that  clearly  demonstrate  such  quantum  effects  in  the  SL  transport. 
To  resolve  this  puzzle  we  develop  a quantum  transport  theory  that  treats  simultaneously 
electro-phonon  resonances  as  well  as  scattering  induced  lifetime  broadening. 

There  is  a fundamental  problem  that  arises  when  the  collisional  broadening  is  taken  into 
account  within  a consequent  quantum-mechanical  approach.  As  the  eigenenergies  of  the 
system  are  no  longer  sharp,  the  electron  distribution  function  depends  explicitly  on  a time 
variable  even  for  stationary  transport  problems.  This  expresses  the  non-Markovian  charac- 
ter of  the  transport.  For  a non-degenerate  electron  gas  a solution  of  the  Dyson  equation  can 
be  searched  for  by  employing  the  following  extended  generalized  Kadanoff-Baym  (KB) 
ansatz 

G<(  fc,  t)  = -(?>(  fc,  t)f  ( k - \t  |,  f)  , (1) 

where  we  introduced  the  notation 


G^(k,t)  = G^ 


e E 

2ft 


t,t  ) , with  G<(  fc,  t)*  — — G<(  fc,  —t). 


(2) 


G ~ are  the  closed-time  path  Green’s  functions.  If  the  collisional  broadening  plays  only  a 
minor  role,  the  explicit  time  dependence  in  the  distribution  function  disappears  (/(  fc,  t)  — ► 
/(  fc))  and  Eq.  (1)  becomes  identical  with  the  well  known  generalized  KB  ansatz  | ],  which 
has  been  used  in  the  literature  to  study  quantum  transport  in  semiconductors.  Unlike  the 
early  KB  ansatz,  which  has  fundamental  limitations,  the  generalized  KB  ansatz  is  fully 
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consistent  with  the  dynamical  structure  of  the  theory  and  agrees  exactly  with  results  derived 
from  the  Liouville  equation  for  the  density  matrix.  This  ansatz  takes  the  causality  for  the 
time  evolution  of  the  particle  propagator  properly  into  account  and  follows  unambiguously 
from  the  symmetry  properties  of  the  electron  system  in  an  external  electric  field.  If,  however, 
lifetime  broadening  becomes  important,  the  generalized  KB  ansatz  no  longer  solves  the 
kinetic  equation  and  one  has  to  determine  a distribution  function  /$(  fc  — e E \ t \/2Ti,  t ) 
that  depends  explicitly  on  a time  variable.  This  leads  to  additional  complications  because 
a closed  equation  cannot  be  derived  for  the  distribution  function  /$  ( fc,  t = 0)  which  is 
used  to  calculate  the  current  density. 

In  this  paper  we  will  study  the  damping  of  electron-phonon  resonances  in  the  SL  transport. 
We  will  show  below  that  such  resonances  survive  only  when  the  lifetime  broadening  is 
extremely  small.  In  this  case  /(  fc,  t)  is  nearly  independent  of  t and  the  kinetic  equation 
simplifies  accordingly. 

2 Numerical  results  and  discussion 

We  consider  a simple  tight-binding  energy  band  of  the  SL 

e(  fc)  = e(  fcj_)  +£(fcz)  = \ k}  + 4(1  - cos k,d),  (3) 

2m*  2 

where  m*  is  the  effective  mass  and  A the  miniband  width.  G>  is  calculated  from  the  Dyson 
equation  under  the  condition  G<  ~ 0,  which  holds  true  for  a non-degenerate  electron  gas. 
In  quasi-classical  approximation  and  for  elastic  scattering  on  impurities  we  get 
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where  the  damping  function  ,v  (£) 
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is  proportional  to  the  scattering  strength  u2 . K is  the  complete  elliptic  integral  of  the  first 
kind.  The  current  density  is  calculated  from  the  stationary  electron  distribution  function, 
which  is  the  solution  of  a kinetic  equation.  Integrating  by  parts,  we  obtain 
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9/(  fc,  t = 0) 
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(6) 


where  V is  the  volume  of  the  crystal.  The  collision  integral  can  be  introduced  in  Eq.  (6) 
when  9/(  fc,  0) / 9/c-  is  replaced  by  the  right-hand  side  of  the  kinetic  equation.  In  the  case 
of  high  electric  fields  and  weak  scattering  (Qt  > 1)  we  get 
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Fig.  1.  (a)  Field  dependence  of  the  dimensionless  current  density  jz/jzo  with  jz o = 

em*nsa>2r /Inftr'd  for  A/Hcoq  = 1 and  Hcoq/IcbT  = 5.  The  scattering  strength  parameter 
m*u2 /n2  Ad  is  given  by  0.005  and  0.05  for  the  solid  and  dashed  line,  respectively. 

(b)  The  same  as  in  (a)  for  A/fkoo  = 0.5. 

where  the  lateral  electron  distribution  function  has  been  replaced  by  the  Boltzmann  dis- 
tribution. ns  is  the  sheet  density,  F the  coupling  constant,  ojq  the  phonon  frequency, 
= Hwo/ksT,  and  s(s,  s')  = s(e)  + s (£■')■ 

Within  our  approach  we  obtained  a broadened  Lorentzian  energy  conservation  in  Eq.  (7). 
This  approximation  has  the  defect  that  higher  and  higher  energy  states  become  populated 
because  the  Lorentz  curve  falls  off  only  gradually.  To  avoid  this  run-away  effect  the  ex- 
plicit time  dependence  of  the  distribution  function  has  to  be  retained.  In  our  present  analytic 
study,  however,  we  will  not  address  the  details  of  such  an  analysis,  but  present  some  nu- 
merical results  that  reveal  already  the  main  features  of  ICFEs  and  collisional  broadening 
in  the  SL  transport. 

Numerical  results  calculated  from  Eq.  (7)  are  shown  in  Figs.  1 (a)  and  (b)  for  A/hcoo  — 
1 and  0.5,  respectively.  Vertical  lines  mark  the  positions  of  electro-phonon  resonances 
at  E — hwo/led  with  / = 1, 2,  3 (for  d we  used  10  nm).  The  solid  lines  have  been 
calculated  for  the  case  when  the  impurity  scattering  strength  is  much  smaller  than  the 
miniband  width  (m*u2 /tt2  Ad  — 0.005).  In  this  case  weak  current  oscillations  appear, 
which,  however,  are  rapidly  smeared  out,  when  the  impurity  strength  becomes  slightly 
larger.  This  is  shown  by  the  dashed  lines,  which  have  been  calculated  for  m*u2/ji2Ad  — 
0.05.  The  lifetime  broadening  effect,  obtained  from  the  microscopic  model,  seems  to  be 
larger  than  phenomenological  estimates  suggest.  This  is  due  to  the  fact  that  in  Eq.  (7)  both 
energy  integrals  are  affected  by  the  damping  function  s(s). 

Our  numerical  results  demonstrate  that  current  oscillations,  which  are  due  to  ICFEs,  occur 
only,  when  the  lifetime  broadening  is  extremely  small.  This  might  be  the  reason,  why 
quantum  mechanical  current  oscillations  have  not  been  observed  in  high-held  transport 
measurements  in  SLs  until  now. 

In  conclusion  we  believe  that  our  quantum-mechanical  approach,  which  accounts  for 
both  ICFEs  and  collision  broadening,  can  be  used  to  study  the  stationary  transport  in  other 
nanostructure  devices,  too. 
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